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ABSTRACT 

Aims. We study the production of VHE emission in blazars as a superposition of a steady component from a baryonic jet and a time- 
dependent contribution from an inner e~e + beam launched by the black hole. 

Methods. Both primary relativistic electrons and protons are injected in the jet, and the particle distributions along it are found by 
solving a one-dimensional transport equation that accounts for convection and cooling. The short-timescale variability of the emission 
is explained by local pair injections in turbulent regions of the inner beam. 

Results. For illustration, we apply the model to the case of PKS 2155-304, reproducing a quiescent state of emission with inverse 
Compton and synchrotron radiation from primary electrons, as well as proton-proton interactions in the jet. The latter also yield 
an accompanying neutrino flux that could be observed with a new generation km-scale detector in the northern hemisphere such as 
KM3NeT. 
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1. Introduction 

Blazars are the AGNs in which the jet points mainly in the di- 
rection of the line of sight. They exhibit the most extreme high- 
energy phenomena of all AGNs. Their spectral energy distribu- 
tions (SEDs) are characterized by nonthermal continuum spectra 
with a broad low-frequency component from X-rays to y-rays. 
Blazars show rapid variability across the entire electromagnetic 
spectrum. Variability at high energies on timescales of a few 
minutes has been observed for some of them, such as PKS 2155- 
304 (e.g. Aharonian et al. 2006). 

In this work we present a two-component jet model with both 
relativistic leptons and hadrons to explain the high-energy emis- 
sion from these objects. The basic scenario consists of a steady 
baryonic jet launched by the accretion disk, and an e + <?~ beam 
launched by the black hole ergosphere. The quiescent compo- 
nent of the signal is assumed to be produced by the jet, while 
the variable component is due to shocked regions in the inner 
<? + e~ beam. Inhomogeneities and turbulence can be generated 
by Kelvin-Hemholtz instabilities. In Sect. 2 we describe the ba- 
sics of the model. Its application to PKS 2155-304 is presented 
in Sects. 3 and 4 for the quiescent and variable emission, respec- 
tively. In Sect. 5 we focus on the neutrino output expected for 
the same blazar, analyzing the detectability with a next genera- 
tion neutrino telescope such as KM3NeT. We finish in Sect. 6 
with a discussion. 

2. Description of the model 

We assume that matter is captured by the central black hole 
through a dissipationless accretion disk (Kelner & Bogolobov 
2010) and that a fraction of this accreted material is expelled by 
the accretion disk in two oppositely directed jets. An inner beam 



of relativistic electrons and positrons is launched by the spin- 
ning black hole. This two-component setup is similar to the ones 
implemented in several previous models (e.g. Sol et al. 1989; 
Romero 1995; Ghisellini et al. 2005; Boutelier et al. 2008). A 
sketch of the basic elements of the scenario is depicted in Fig.Q] 
The introduction of a model with two components is moti- 
vated by observations as a means to reconciling available data 
with the unified AGN paradigm (e.g. Urry & Padovani 1995), as 
noted by Chiaberge et al. (2000). One-component jet models re- 
quire high bulk Lorentz factors (10-20) when applied to blazars 
and lower values when applied to FR radiogalaxies. Including of 
a second component serves to solve this discrepancy, for exam- 
ple, through the combination of a fast spine and a slower (but still 
relativistic) layer, so that the emission at small viewing angles is 
dominated by the fast component, whereas the slow component 
dominates for larger angles. 

Some additional observational facts also favor the two- 
component interpretation. For instance, the presence of a fast 
spine surrounded by a slower outflow can be inferred from the 
observed limb-brightened radio morphology of the radiogalax- 
ies jets (Giroletti et al. 2004). On the other hand, whereas VLBI 
observations of powerful TeV BL Lacs suggest that the pc-scale 
jets move slowly (Piner & Edwards 2004; Giroletti et al. 2004), 
rapid variability of the TeV emission implies that, in the region 
where this emission is produced, the jet should be extremely rel- 
ativistic (e.g., Dondi & Ghisellini 1995; Ghisellini et al. 2002; 
Konopelko et al. 2003). 

On the theoretical side, the presence of both a spinning 
black hole and an accretion disk would unavoidably yield the 
launching of outflows through the BZ process (Blandford & 
Znajek 1977) and the BP process (Blandford & Payne 1982). 
The coexistence of both mechanisms has been investigated by 
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Meier (2003), among others. Whereas the BZ process yields the 
launching of a Poynting flux that gives rise to a leptonic beam, 
the BP process can generate a baryon-rich jet launched centrifu- 
gally from the inner accretion disk (e.g. Komissarov et al. 2007; 
Sadowsky & Sikora 2010). 

The stability and structure of this type of jets have been dis- 
cussed, for instance, in Hardee (2007), Narayan et al. (2009), 
and Perucho (2012). For general discussion of fluid instabilies 
see, e.g., Shore (2008). 




Fig. 1. Basic elements of the model. See the text for details. 

For the jet, we adopt the model developed by Reynoso et al. 
(2011). Equipartition between jet kinetic energy and magnetic 
energy takes place at zq = 50R g from the black hole, the jet 
half-opening angle is £j, the jet Lorentz factor is Tj at zq, and the 
viewing angle is l s . The bulk kinetic power of the jet at zo is a 
fraction of the Eddington power, 



(kin) 



<7j£-Edd = (Fj - l)m s c 



and the cold particle density in the jet is 



'J - 2 2' 



(1) 



(2) 



with^j = tan^j. The magnetic field varies along the jet according 
to 

B(z) = Bo(— J withm e (1,2), 

and the Tj increases slowly as magnetic energy density decreases 
along the jet. 

The e + e~-beam, with a half-opening angle < £j, is con- 
fined by the jet by requiring that the pressure is less than in the 
jet: P\, — ilbPj, where t]\, < 1. The initial Lorentz factor of the 
beam is F b 0) at distance z b °' < zo from the black hole, where 
equipartition with the beam magnetic energy holds. This condi- 



tion is used to fix B 



.(()) 



For z > z b 0) in the beam, we assume 



Bb - B\ 



.«)) 



with nit, 6 (1,2), 



(3) 



so that the Lorentz factor of the beam also increases gradually 
along it. This behavior of the magnetic field also enables Kelvin- 
Hemlholtz instabilities to develop, as we discuss below in Sect. 



|U where we deal with the production of variable emission in the 
beam. As mentioned, this contribution is added to the one from 
the quiescent jet, which is discussed in next. 

3. The quiescence emission from the jet 

The quiescent state of the emission can be reproduced using a 
lepto-hadronic model for the jet. Acceleration of relativistic elec- 
trons and protons takes place at z acc > zo, where the power in- 
jected in these primary particles is L e and L p = aL e , with a > 0. 
The total power in the injected, L re \ = L e ( \ + a), is a fraction of 
the bulk kinetic power of the jet, L re i = q , re iLj kin ' > (see Reynoso et 
al. 2011). The distribution of particles in the steady jet, N(E,z) 
in units GeV^'cirT 3 , is governed by an inhomogeneous transport 
equation with cooling and convection: 



d(Tj Vi N(E,z)) + d(b(E,z)N(E,z)) + N(E,z) 



dz 



BE 



Tdec(E) 



Q(E,z), (4) 



and 



where z is the distance to the black hole in the AGN frame, 
b(E,z) = -§. 

This equation is solved in the jet co-moving frame, and the 
injection or source term is taken as 



Q(E,z) = K, 



- i Zacc \^ 



E~ s exp 



-.(max) 



(5) 



where the cut-off energy £^ max) is determined by equating the 
acceleration rate to the cooling one. The particle injection is nor- 
malized in the AGN frame, that where the central black hole is 
at rest, as 



Li = f dE' \ dQf \ dV E 

JAE J4n JAV 



, d Qj(E\z) 
dCl' ' 



where 



and 



dQ£E\z) 
dD,' 



E' 2 - m 2 c 4 



E 2 - m 2 c 4 



1/2 



Qi(E,z) 
An ' 



E = r^E" - Bj cos 9 yj E' 2 -rrvrc 4 ^. 



The acceleration rate is given by 



t;l(E,z) 



ceB(z) 



(6) 



(7) 



(8) 



(9) 



with an efficiency rj, and the loss rate corresponds to synchrotron 
emission, inverse Compton interactions, pp and py interactions, 
and adiabatic cooling (see Reynoso et al. 2011). In Fig. [2] we 
show the acceleration and cooling rates for electrons and protons 
at the position z acc in the jet. The list of assumed and derived 
parameters for the jet is shown in Table 1 . The values for the 
black hole mass and redshift are taken as Mbh = 1O 9 M , and 
z rs = 0.1 16, respectively, following Aharonian et al. (2007), and 
the column density of HI is taken as Nh = 1.3 X 10 20 crrT 2 after 
Lockman & Savage (1995). The rest of the parameters in Table 1 
corresponds to the above described jet model, which is discussed 
in detail in Reynoso et al. (201 1). These parameters are fixed to 
obtain a SED that fits the quiescence emission represented by the 
same set of observational data used in Aharonian et al. (2009). 

The SED is dominated by synchrotron and inverse Compton 
(IC) emission of the electrons, and pp orginated VHE gamma- 
rays, as shown in Fig. [3] In our calculations, the VHE energy 
emission has been corrected by the effect of EBL absorption fol- 
lowing Dominguez et al. (2010). 
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Fig. 2. Acceleration and cooling rates for primary electrons and protons in the jet at z : 



Table 1. Model parameters for the quiecence emission from the 
jet of PKS 2155-304 



PKS 2155-304, quiescence SED 



Parameter 


Description 


Value 


M bh 


black hole mass 


10 y M o 


Is 


ratio 2Zi km) /£Edd 


0.1 


rj(zo) 


bulk Lorentz factor of the jet at zo 


10 


h 


viewing angle 


1° 




jet's half-opening angle 


1.5° 


<?rel 


jet's content of relativistic particles 


0.1 


a 


hadron-to-lepton power ratio 


50 


Zo 


jet's launching point 


50 R gl 


q m 


magnetic to kinetic energy ratio at z acc 


0.02 " 


Zacc 


injection point 


575 R g 


A Zj 


size of injection zone 


z acc tanfj 


m 


index for magnetic field 


1.8 


s 


spectral index injection 


2.1 




acceleration efficiency 


10-" 


ri(min) 


minimum proton energy 


2GeV 




minimum electron energy 


0.9 GeV 


N H 


column density of HI 


1.3 x 10 20 crrT 2 


Zis 


redshift 


0.116 



I 1 I 1 I 1 I 1 I 1 I 1 I 1 I 1 I 1 

+ Data (quiescent state) — — - IC 

^— ^— e-syn — ■— pp (w/EBL abs) 

^— ^— p-syn pp (no EBL abs) 

— — — py (no EBL abs) 




10 4 10 6 10 s 10 1 ' 
E v [eV] 



R„ = 1.48 x 10 14 cm 



Fig. 3. SED corresponding to the quiescent state of PKS 2155- 
304. The contributions shown are electron synchrotron (black 
solid line), IC (dashed line), pp (dashed-dotted line), proton syn- 
chrotron (gray solid line), and py (gray dashed line). 



4. Variable emission from the beam 

The variable emission is assumed here to have its origin in 
shocked regions in the internal e + <?~ beam that can arise, e.g., 
due to Kelving-Helmholtz (KH) instabilities. These instabilities 
will be able to develop in the beam if the magnetic field is below 
a critical value (e.g. Romero 1995): 

jAnn h m e c 2 {Tl - 1) 
B h < Bf 1 = J . (10) 

I b 

In Fig.|4]we show the beam magnetic field and the critical value 
B™, considering that r£ = 5 at = 5R g . It can be seen from 
this plot that for z > 10 16 cm KH instabilities start to develop 
and create inhomogeneities if the pinching modes dominate (see 
Romero 1995). The rarefactions produced in the beam by the in- 
stabilities act as obstacles for the fast plasma, and strong shocks 
are expected to appear at z ~ 10 cm. These shocks reacceler- 
ate the pairs (Araudo et al. 2010) producing a local injection of 
relativistic particles. 
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Fig. 4. Magnetic field in the beam (red line) and critical magnetic 
field below which KH instabilities can develop. 
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A shock front originating in the beam at a position z; with 
a Lorentz factor F s will form a post-shocked region of a certain 
size 



ZiXb 



with^b = tan^t,, and particles are injected during a time 

Az 



At 



inj, obs 



(ID 



(12) 



as seen in the observer frame. 

The distribution of the injected pairs is computed in the post- 
shock frame, where the population is assumed to be isotropic. 
The corresponding source term is a power law with an exponen- 
tial cutoff at the maximum Lorentz factor y mdx : 



q e (y, t) = H(t - t on )H(t oB -t)q y Ib exp 



y 



Tir 



(13) 



is where H is the Heavyside step function, t on the time when the 
injection is switched on, and the time when it is switched off is 



At 



to« - ton + 



inj, obs 



r ps (l + Zrs) 



(14) 



both corresponding to the post-shock frame. We note that, for 
example, the onset time is transformed to the observer frame as 



^on, obs — l~p 8 (l Zrs)^on- 



(15) 



The constant go is found through normalization taking into 
account the particle density in the post-shock frame: 



X'oir r"*> 
dt J dyq e (y,t). 



(16) 



This can be related to the beam density n by the compresion 
factor £ as 

» ps = <r"b, with { = ( yT / s _ + ^ , (i7) 

where y = 4/3 is the polytropic index, 

r; s = r ps r b (i-/? ps /j b ) 

is the Lorentz factor of the shock in the undisturbed beam frame, 
and the Lorentz factor of the post-shock (ps) region is y ps * 
y s / V2. The density of the beam, in turn, is fixed as a fraction A 
of the jet density of Eq. [2] 



(18) 



The maximum Lorentz factor of the locally accelerated pairs 
Tmax is obtained from the balance between the acceleration rate 
and the cooling rate, where the latter is due to synchrotron and 
synchrotron self-Compton (SSC) interactions. 

We solve the following kinetic equation to obtain the distri- 
bution of pairs as a function of time n e (y, t) : 



dn e d 
where the loss term is given by 

Itlloss = Itlsyn + ItlsST • 



(19) 



(20) 



An iterative process is implemented to work out n e (y, t) from 
this equation. As a zero-order approximation, we adopt the so- 
lution corresponding to no SSC interactions, just synchrotron 



losses. Then, the solution is succesively improved by comput- 
ing the SSC cooling term using the previus solutions found for 
n e (y, t). This cooling rate is computed in the Thompson regime 
to save computing time, since this simplification was tested not 
to affect the final result. The self-synchrotron Thomson (SST) 
cooling rate is then (e.g. Schlickeiser & Lerche 2007): 



X 



where 



ItlsST = A y dyy n e (y, t), 



(21) 



Ao = 3 CCTT 8^ 
The synchrotron cooling rate is 

Irlsyn = - y ■ 



6nm e c 



(22) 



The acceleration rate is worked out using Eq. (O along with the 
magnetic field corresponding to the post-shock frame, 



B ps (Zi) = B b ( Zi ) ^F+l. 



(23) 



For any injection point z, in the beam, we solve the kinetic 
equation, and work out the resulting synchrotron and SSC emis- 
sion. In the local AGN frame, the emissivity is given by 



dQ' y {E',t') _ dQ y {E,t) 
dQ.' ~ dQ. ' 



(24) 



where the Doppler factor is D = r~ (1 - /? ps cos , E is the 
photon energy in the ps frame, and E' — D E is the photon 
energy in the local AGN frame. In the case of synchrotron radi- 
ation, the emissivity in the ps-frame is 

dQr,syu(E,t) 1 ( l- e -™<M \ f° 

~^T- - M W *,o J I, * p ^ • f) ' (25) 

where tssa is the optical depth corresponding to synchrotron- 
self absorption, and the synchrotron power per unit energy emit- 
ted by the electrons is given by (Blumenthal & Gould 1970): 



fsy„(£,7^) : 



V2e 3 g ps E 
m e c 2 h E cr J E/Ea 



d{K 5n (0, 



(26) 



where X5/3 is the modified Bessel function of order 5/3 and 



p - 



V6/ieB ps 



-T 



Anm e c 

The soft photon density in the ps frame is 

and the IC emissivity in the this frame is 
r 2 c C 

Q y ,ic(E,t)=^- J e 



(27) 



rL n vh (E ph ,t) 
dE p h 



j P h 

X r a dy'^Apiq). (28) 



r 



Here r e is the classical electron radius, and we integrate in the 
target photon energy £p n and in the electron Lorentz factors y' 
between 



2m„c 2 



E E 2 

1 7 and y„ 

4£ ph Amy Yn 



m e c 2 (l-^Y 
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The function F(q) is given by 



Table 2. Beam parameters 



F{q) = 2q]nq + (l + 2q){l - q) + 1(1 - qij^, (29) 
with T e = AE v hy' l(m e c 2 ) and 



1 = 



r e £ , p h(l - E/E p h) 



The differential photon flux is first obtained in the AGN local 
frame as 



, r -i 
7 exp[-r yr J, (30) 



— = (ttz 2 tan 2 fl,Az) 

d£' ' % dQ.' dA 

where r yy = rS + ?w is the yy optical depth due to the jet and 
the beam photons. The latter is estimated as 

t™(E', t') = ( Zi tan^ b ) j" dx(l - x) 

x I ^ ^phMph^ph, t)<r 7 y(E, E pb ), (31) 

and r yr is obtained from a similar expression. 
If dc is the comoving distance, then 

da' _ i 
Ja> ~ 4' 

and taking into account that the observed energy and time are 
transformed as 



E y = 



E' 



1 +Zrs 



and f obs = f'(l +z rs ), 



the differential photon flux arriving on Earth from a source at 
redshift Zt$ is 

d%(E y , f obs ) (1 + Zrs ) 2 d% [e 7 (\ + Zrs ), ^-) 



dE y 



d[ 



dE' 



exp-^EBL^y.Zrs)], (32) 



where d^ is the luminosity distance and tebl is the optical depth 
due to yy absorption on the extragalactic background light. In 
these terms, the usual amount vF v in units erg ctrr 2 s~' is ex- 
pressed as Ey'^r-. A light curve can be obtained simply by inte- 
grating Eq. ( l32l on the desired energy range for different times. 
In Fig. [5]we show the light curves that can be obtained if the in- 
jection point in the beam is Zi = 2.3 x 10 16 cm for different shock 
Lorentz factors T s and sizes of the injection zone Az — R h /b. The 
values adopted for beam parameters are shown in Table [4] More 
complex light curves can be obtained, for instance, by succesive 
shock injections at the same position with T s , b, and f on chosen 
randomly in the ranges: 



T s e {17,18,19} 
b e {50,60,70,80,90} 
t on e [1..15]min. 



(33) 
(34) 
(35) 



Some example light curves obtained in this way are plotted in 
Fig.|6]with the parameters shown in Table[3] 

In particular, in view of the flaring activity detected in 2006 
(Aharonian et al. 2007), we see that a simillar lightcurve can be 



Parameter 

-(b) 
p(0) 
1 b 
A 

& 



Description 



launching point [R g ] 
Lorentz factor at 
beam to jet density ratio 
half-opening angle [°] 



Value 



5 
5 

0.1 
0.5 



Table 3. Randomly chosen values for the parameters corre- 
sponding to Fig|6] 



Label 


r s 


b 


f on [min] 


al 


17, 18, 19, 19, 18 


80, 80, 60, 50, 80 


5.3,8.1,5.5,11 


al 


17, 17, 18, 18, 19 


50,90,50,50,50 


6.3,3.6,11,12.4 


a3 


18,17,19,19, 17 


90, 80, 80, 80, 80 


4.2, 2.6, 12.3,9.4 


a4 


17, 18, 17, 18, 19 


90, 70, 80, 90, 70 


6,7.2, 1.7, 12.3 


a5 


17, 18, 19, 18, 18 


90, 60, 80, 60,90 


3, 10.7,7, 10.1 


a6 


18, 18, 17, 18, 17 


60, 60,70,90, 80 


9.1, 12.9,8.1,2.2 


al 


17, 17, 17, 18, 18 


80, 70, 80, 50, 60 


11.6, 2.8,4.4,2.1 


a8 


17, 19, 17, 18, 18 


80,70, 80,50, 60 


7.6,5.2,3.5,9.6 


a9 


19, 17, 17, 19, 17 


50, 60, 60,70, 60 


11.7,4.7,5.5,7.6 



PKS 2155-304 (2006 flare 




t [min] 



Fig. 7. Lightcurve of gamma-rays above 200 GeV (blue line), 
generated as the superposition of six injections in the beam 
(lei Jc2,lc3,lc4,lc5,lc6, in gray-dashed lines), as compared to 
data of the 2006 flare of PKS 2155-304. 



obtained for the flux of gamma-rays above 200 GeV by adding 
up six different injections at a fixed position in the beam, Zi = 4x 
10 16 cm, each being switched on at apropriate times f on , as seen in 
Fig-El The values of the different parameters are shown in Table 
[4] The multiple shock scenario resembles the one proposed to 
explain the rapid variability in gamma-ray bursts (GRBs), e.g., 
Kobayashi et al. (1997). The whole SED evolution shown in Fig. 
8 shows that multiwavelength observations, including studies of 
rapid variability in radio, X-ray, and y-ray bands, can be used to 
test the proposed model. 

In Fig. [8] we show for illustration two instantaneus contri- 
butions to the SED due to the shocked beam at t\ — 40min and 
?2 — 50min. The acceleration and cooling rates corresponding 
to 1 1 are shown in Fig. [9] where it can be seen that the Lorentz 
factors of the electrons are ~ 10 6 in the post-shock frame. The 
inverse Compton interactions with target photons created in the 
jet was verified to be unimportant for the adopted value of zu 
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Fig. 5. Light curves of gamma-rays above 200 GeV produced by e + e injection in the beam at zi = 2.3 x 10 cm with r s = 
{14, 16, 18, 20} in the left panel, and with T s = 20 and b = {50, 70, 100} in the right panel. 




„ 5xl0 9 
fc 4xl0" 9 



<x» 

o 2xl0 9 

fN 

A 

uT lO" 9 





1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 




a7 




1 N 


I 














A / « i" 


i 






r \ 1 '» <! 1 








t" \ 1 i \ \i 








F \ 1 i m\ 








r 1 1 1 n 1 1 1 1 1 1 1 1 1 1 ih+4. 


444 1 II 


-J 




10 20 30 40 50 60 
t bs [min] 



10 20 30 40 50 60 
tobs [min] 



10 20 30 40 50 60 
tobs [min] 



Fig. 6. Light curves of gamma-rays above 200 GeV produced by e + e injection in the beam at Zi = 2.3 x 10 16 cm with the randomly 
chosen values for T s , b, and t oa of Table|3] 



5. Neutrino emission 

The acceleration of protons in the jet and their interactions with 
synchrotron photons and cold jet protons lead to the production 
of VHE gamma-rays after ir decays, as seen in the SED of Fig. 
[3] These py and pp interactions also give rise to charged pions 
that decay to muons and neutrinos. The stationary distributions 



of pions and muons in the jet, N n (E,z) and Nf,(E,z), are found 
using the transport equation Eq. [4] The corresponding neutrino 
emissivity from direct pion decay can be computed in the jet 
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Table 4. Model parameters for the variable emission from the beam of PKS 2155-304 



Parameter 


description 


lcl value 


lc2 value 


lc3 value 


lc4 value 


lc5 value 


lc6 value 


r s 


shock Lorentz factor 


13 


10.5 


13 


13 


11.5 


11.5 


% 


acceleration efficiency 


io- 2 


IO 2 


io- 2 


io- 2 


io- 2 


io- 2 




power law injection index 


2.1 


2.1 


2.1 


2.1 


2.1 


2.1 


Zi 


position of injection at the beam [cm] 


4x 10 16 


4x 10 16 


4x 10 16 


4x 10 16 


4x 10 16 


4x 10 16 


Az 


size of injection zone 


z^b/50 


z<*b/50 
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Fig. 8. Instantaneus SEDs originating in the beam at t\ ^ 40min, 
and ?2 - 50min corresponding to the plot of Fig(7] The quiescent 
SED is also shown (gray dashed lines). 
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Fig. 9. Acceleration and cooling rates in the post-shock frame 
corresponding to t \ ^ 40min in the lightcurve of FigJT] 



comoving frame following Lipari et al. (2007): 



<?.-„<£.:) = f dE n T-^E n -)N n (E n ,z) 



x *_^i), (36) 

E„(l - r„) 



with x = E/E n and = 2.6 x 10 s. The contribution from 
muon decays {[T — > e^v e v^, — > e + v e y tl ) can be calculated as 

4 poo 
,= 1 J E ^ 



5 , 2 4 3 

3x + -x 

3 3 



(37) 



In this expression, x = E/E^, \i\ 2 - H L ' + , T^a = 2.2 x 10 6 s, 

andyu 3 ,4 = H R ' + - 

In the local AGN frame, the differential neutrino emissivity, 
in units (GeV^'cm^sr^s" 1 ) is given by 

dQ' v (E',z) = Q^(E,z) + Q^(E,z) 

d& An ' 1 ' 

where z is the distance from the black hole to the position in the 
jet in the AGN frame and D - Ft 1 ( 1 cos ij) -1 . The differential 
neutrino flux is then 

r ,dQ' v (E',z') 



dE> dl J z -> 



dV 



dQ.' 



(39) 



and in the observer frame, since the redshifted neutrino energy 
hE v = E'/(l+Zn), 



d® v (E v ) d%(E v (l+Zrs)) 



dE v 



dE' 



(40) 



On Earth, neutrino telescopes are used to observe this flux in 
the presence of the atmospheric neutrino background, ^nr§$ , 
which introduces more events depending on the size of the an- 
gular search bin of the detector, 



Afihin = 27T 1 - COS 



AObh 



This depends on the resolution of the detector Aflbin, which is 
expected to be below 1° for KM3NeT. In the case of PKS 2155- 
304, we compare in Fig[l0]the predicted neutrino flux with the 
atmospheric background for A#bj n = {0.5°, 1°). 

The signal to be searched for in neutrino telescopes such as 
ANTAPvES or KM3NeT, can be calculated as 



N S (E V > £t mm) ) = T, 



obs I 

J Eg 



<min) dEy 



(41) 



where Af{E Y ) is the neutrino effective area of the detector and 
r o bs is the observation time. These signal events are to be de- 
tected among the background events caused by the flux of atmo- 
spheric neutrinos, which are given by 



N b (E v > E v mm >) = TobsAQi 



bin 



Xmi„: 



dE v dQ. 



(42) 
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Fig. 10. Predicted muon neutrino differential flux weighted by 
the squared neutrino energy as compared to the flux of atmo- 
spheric neutrinos for A#bj n = {0.5°, 1°). 
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Fig. 11. Muon neutrino events from PKS 2155-304 and atmo- 
spheric background for A#bj n = {0.5°, 1°) as a function of the 
observation time. 

Considering the neutrino effective area expected for this de- 
tector, we compute the predicted number of signal events of 
(Vft + v^) from PKS 2155-304 as a function of the observation 
time, and compare it with the background events corresponiding 
to Af?bin = 1° and 0.5°. The result is shown in Fig.fTTIand corre- 
sponds to neutrinos above = 1 TeV, yielding many events 
in a few years. 

We can proceed to compute the cumulative probability that at 
least one of the total events corresponds to the signal, and not to 
the background, as a function of the observation time. This can 
be done by assuming that the ocurrence of the signal event and 
backgroungd events N s and N b follow Poisson distributions. As 
a result, the probability that out of a total of events N tot = N s +Nb 
we have at most N\, < A tot can be expressed as 

N ^ A*exp(-A b ) 
P(N b < AW) = X - ± —r, ■ (43) 

We show this cumulative probability in FigQ~2]for both Af9bm = 
{0.5°, 1°), where it can be seen that if PKS 2155-304 actually 
emits neutrinos at the level suggested here, a neutrino telescope 
like the planned KM3NeT would be able to ascribe at least one 
event to this blazar with a probability very close to one after an 
observation period of four years. We point out here that a reso- 
lution of 0.5° is actually conservative (e.g. Kappes et al. 2007, 
Bersani 2012). 



0.6 1 ' 1 ' 1 1 1 ' ' ' 1 1 1 ' 1 ' ' 1 1 

2 4 6 8 10 

T[yr] 

Fig. 12. Cumulative probability that at least one event corre- 
sponds to the neutrino signal from PKS 2155-304 as a function 
of the observation time. 

6. Discussion 

We have implemented a two-component model to study the 
emission of the blazar PKS 2155-304. The quiescent state of 
electromagnetic emission was associated with the contribution 
produced in a heavy barionic jet, in which both electrons and 
protons can be accelerated. Specifically, this low state of emis- 
sion is dominated by synchrotron and inverse Compton inter- 
actions of the primary electrons, and also by pp interactions in 
the jet. The time-dependent contribution is produced by multiple 
shocks in an internal electron-positron beam, surrounded by the 
jet. We showed that if shocks are injected in this beam, as would 
be expected as a consequence of Kelvin-Hemlholtz instabilities, 
then a variable emission can be generated via SSC interactions, 
giving rise to gamma-ray light curves similar to those observed 
by HESS. 

Other important output predicted by the present model is in 
the form of VHE neutrinos produced by pp interactions. If this 
accompanying flux remains constantly produced over a period of 
four years, the detection of at least one event from PKS 2155-304 
should be guaranteed in a detector like KM3NeT with a proba- 
bility very close to one. This would be a very important piece 
of evidence for hadronic acceleration in the source. Conversly, 
if the neutrinos are not detected at this level, this would be evi- 
dence of a leptonic-dominated source. 
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